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Models that lead to a cosmological stiff fluid component, with a density ps that scales as a -6 , 
where a is the scale factor, have been proposed recently in a variety of contexts. We calculate 
numerically the effect of such a stiff fluid on the primordial element abundances. Because the 
stiff fluid energy density decreases with the scale factor more rapidly than radiation, it produces a 
relatively larger change in the primordial helium-4 abundance than in the other element abundances, 
relative to the changes produced by an additional radiation component. We show that the helium-4 
abundance varies linearly with the density of the stiff fluid at a fixed fiducial temperature. Taking 
psio and pmo to be the stiff fluid energy density and the standard density in relativistic particles, 
respectively, at T — 10 MeV, we find that the change in the primordial helium abundance is well-fit 
by AY P — 0.00024(psio/pflio). The changes in the helium-4 abundance produced by additional 
radiation or by a stiff fluid are identical when these two components have equal density at a "pivot 
temperature", T», where we find T» = 0.55 MeV. Current estimates of the primordial 4 He abundance 
give the constraint on a stiff fluid energy density of psio/pmo < 30. 



I. INTRODUCTION 

In recent years, a host of cosmological observations 
have provided an increasingly precise picture of the con- 
stituents of the universe. The baryon density has long 
been known to provide roughly 5% of the critical density; 
earlier estimates from Big Bang nucleosynthesis [l| have 
been confirmed by CMB observations from WMAP 0. 
Cosmological data from a wide range of sources including 
type la supernovae [|| , the cosmic microwave background 
_baryon acoustic oscillations [3 , cluster gas fractions 
and gamma ray bursts 0, H| seem to indicate that 
about 72% of the energy density of the Universe is in 
the form of an exotic, negative-pressure component called 
dark energy. In addition, a slew of evidence from mod- 
ern sources, including weak [9j and strong II Oil lensing, 
the bullet cluster [ll| , Large Scale Structure [l2| , as well 
as supernovae and the CMB, have confirmed earlier in- 
dications from rotation curves [l3T - fll| that another 23% 
of the energy density of the Universe is in the form of 
weakly interacting matter called dark matter. (See Ref. 
Hi. for a review of dark energy and [l7| for a review of 
dark matter). 

While the existence of all of these components is rea- 
sonably well-established, the existence of other exotic flu- 
ids is not ruled out by the current data. For example, 
several models predict the existence of "dark radiation" 
either early or late in the history of the Universe (see e.g. 
[l8| for an example of such a model and references therein 
for others). Another exotic fluid which arises in various 
models is a "stiff fluid", i.e., a fluid with an equation of 
state parameter ws = Ps/ps — 1- This is the largest 
value of w consistent with causality, since the speed of 
sound of this fluid equals the speed of light. Such mod- 
els were apparently first studied by Zeldovich 19]. The 
Friedman equation for such a fluid implies that its energy 
density ps varies with the scale factor a as: 

Ps cx a' 6 . (1) 

In recent years, a variety of models have been proposed 



that produce a stiff cosmological fluid: 

a. Kination: A kination field is a scalar field whose 
energy density is dominated by kinetic energy. A period 
of kination can follow a period of inflation, and was first 
studied in the context of electroweak baryogenesis [2(| 
[2lT |. Its impact on reheating as well as on the freeze-out 
of dark matter particles has been studied in Refs. [22l426j j . 
Kination fields have also been applied to unify inflation 
and dark energy by using the same scalar field for both 
inflation and quintessence. In these models inflation ends 
with a period of kination as the inflaton receives a "kick" , 
and the same field later on plays the role of quintessence 

b. Interacting dark matter: In models with a warm 
self-interacting dark matter component, the elastic self- 
interactions between the (scalar boson or fermionic) dark 
matter particles can be characterized by the exchange of 
vector mesons via minimal coupling. For these models 
the self-interaction energy can be shown to behave like a 
stiff fluid [H. 

c. Hofava-Lifshitz cosmologies: Stiff fluids also oc- 
cur in certain cosmological models based on the recently 
proposed Hofava-Lifshitz gravity, a power-counting 
renormalizable and ultraviolet-complete field theoretic 
quantum gravity model based on anisotropic scaling of 
the space and time dimensions [29-32]. In the original 
formulation of this theory a "detailed balance" condition 
was imposed as a convenient simplification 31|. The va- 
lidity and usefulness of the detailed balance condition 
have subsequently been discussed extensively (see e.g . 
[33|,|34j]) as well as the consequences of relaxing it (35l - [38| . 
Stiff fluids arise in models in which the detailed balance 
condition has been relaxed. Cosmological models based 
on Hofava-Lifshitz gravity have been studied extensively 
(see e.g. [39j - l4^ |]) and observational constraints on such 
models, including the stiff-fluid cases, were considered in 
[45M47] ]. However, it must also be noted that the theoret- 
ical foundations of Hofava-Lifshitz gravity are still under 
debate (see e.g. HHsl). 
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d. Non-singular cosmological models: Stiff fluids 
have been also found to show up as exact non- sing ular 
solutions in inhomogeneous cosmological models [57H60J ■ 

Given the recent flurry of interest in such models, it 
is clearly useful to derive precise limits on the density 
of a stiff fluid in the early universe. Because the den- 
sity of a stiff fluid decays more rapidly than either ra- 
diation or matter, the effect of the stiff fluid on the ex- 
pansion rate will be the largest at early times. Thus, 
the strongest limits on the density of such a fluid come 
from Big Bang nucleosynthesis (BBN), which remains the 
earliest cosmological process whose evolution can be de- 
termined with high precision. (Note that some attempts 
have been made to constrain the expansion rate of the 
universe prior to BBN from the relic dark matter abun- 
dance (^mmiUHl. While the freeze-out of the dark 
matter does occur at a much earlier time than BBN, the 
exact model for the dark matter is less precisely deter- 
mined) . 

Previous discussions of stiff fluids have usually quoted 
BBN limits on the expansion rate at a fixed tempera- 
ture (typically T ~ 1 McV) and used these limits to 
constrain the density of the stiff fluid. Here, we numer- 
ically evolve the element abundances in the presence of 
the stiff fluid to derive the exact dependence of the ele- 
ment abundances on the stiff fluid density. The stiff fluid 
investigated here resembles a special case of the mod- 
els examined by Masso and Rota |65[ who investigated 
the effect of an arbitary additional density of the form 
A(T/0.1 MeV) 7 , although their 7 = 6 case is not identi- 
cal to the model considered here, since T does not scale 
as 1/a through the era of e + e~ annihilation. Here, we 
examine the stiff fluid case in more detail, and we ex- 
ploit the fact that WMAP now provides an independent 
determination of the baryon-photon ratio @], effectively 
eliminating a degree of freedom from BBN and allowing 
for better constraints on the stiff fluid. (See also related 
work on a somewhat different variant model in Ref. 66] ) . 
We present our calculations in the next section, and dis- 
cuss our limits in Sec. III. 



II. EFFECT OF A STIFF FLUID ON BBN 

BBN has long been used to constrain additional energy 
density in the early universe (for recent reviews, see Refs. 
[67l - l69j ). The expansion rate H is given by 



neutrons and protons, 



H 2 



8vrG 



(2) 



where p is the total density, so any additional contribu- 
tion to p increases the expansion rate and changes the 
resulting element abundances. 

At high temperatures (T Si 1 MeV) the rates for the 
weak interactions which govern the interconversion of 



n + v e <-> p + e~ , 
n + e + <H> p + v e , 

n o p + e~ + v e 



(3) 



are larger than the expansion rate, H, and the neutron- 
to-proton ratio (n/p) tracks its equilibrium value. As the 
universe expands and cools, the weak interaction rates 
drop below the expansion rate, and n/p freezes out at 
T ~ 1 MeV. Between T ~ 1 MeV and T ~ 0.1 MeV, the 
neutrons undergo free decay, and then at T ~ 0.1 McV, 
nucleosynthesis proceeds to fuse the remaining neutrons 
with the protons to produce heavier elements, primarily 
4 He, but also trace amounts of 2 H, 3 He, and 7 Li. The 
final He abundance is most sensitive to the expansion 
rate near T ~ 1 MeV, when the neutron/proton ratio 
freezes out, while the other element abundances are more 
sensitive to the expansion rate near T ~ 0.1 MeV, when 
fusion into heavier elements begins. 

The dependence of the different element abundances 
on the expansion rate at various temperatures was ex- 
plored quantitatively by Bambi et al. [7(J, who derived 
"response functions" that show the change in the abun- 
dance of each nuclide as a function of a change in the 
expansion rate at a given temperature. [Note that we 
use these response functions only to gain insight into the 
effects of the stiff fluid on the element abundances; our 
actual calculation of the element abundances utilizes a 
full numerical integration of the BBN equations, as dis- 
cussed below]. Consider an additional source of energy 
density which changes the expansion rate, H(T), by an 
amount AH{T), Then Bambi et al. argued that, for a 
given value of the baryon-photon ratio 77, the change in 
a given nuclide abundance, AXi, is given by 



Al, 



Qi(T) 



AH(T) dT 
H(T) ~T ' 



(4) 



where Qi{T) is the response function for a given nuclide, 
derived numerically for several elements of interest in Ref. 
[rol | . As expected, the response functions for deuterium 
and 7 Li are strongly peaked near T ~ 0.1 MeV. In con- 
trast, the response function for 4 He is broadly distributed 
between 1 MeV and 0.1 MeV, with two peaks of roughly 
equal magnitude at these two temperatures 17(1 . The 
first corresponds to n — p freeze-out, and the second to 
the onset of fusion. The latter affects the 4 He abundance 
primarily through the influence of free-neutron decay; the 
earlier fusion begins, the more undecayed neutrons re- 
main to form 4 He. 

Now consider a model, such as the one examined here, 
with some additional source of energy density, Ap(T). 
As long as Ap(T) <C p{T), where p(T) is the standard 
energy density, we can write equation Q as 



AXi 



Qi{T) 



Ap(T) dT 
W) ~T' 



(5) 
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Because of the sensitivity of the element abundances 
to the density at different epochs, the relative changes 
in the different element abundances will be different 
for different functional forms of p(T). However, con- 
sider an arbitrary functional dependence of the form 
p(T) = p(T )f(T/T ), where p(T ) is the density at some 
fixed fiducial temperature To, while /(T/Tq) is an arbi- 
trary function subject only to the constraint /(l) = f. 
Thus, f(T/T ) parametrizes the dependence of p on T, 
while p(T ) fixes the overall amplitude of the density. 
The important point is that as long as Eq. J5]) is a good 
approximation, and once f(T/Tg) is fixed, the abundance 
of each nuclide will vary linearly with p(Tq), regardless 
of the functional form of /(T/Xb). 

This linear dependence is seen in the case of addi- 
tional relativistic energy density scaling as Apu oc a~ 4 . 
Parametrizing this energy density in terms of the number 
of additional two-component neutrinos, AN V , the change 
in the primordial 4 He mass fraction, AY p , is well approx- 
imated by 0| 

AY P = 0.013 AN„. (6) 

Given the different sensitivities of the element abun- 
dances to the expansion rate at different temperatures, 
it is clear that the change in the element abundances pro- 
duced by a stiff fluid will differ from that produced by 
additional relativistic energy density. However, if we con- 
fine our attention to a single element (such as 4 He), then 
we expect the overall abundance to scale linearly with 
the value of the stiff energy density at a fixed fiducial 
temperature. 

We model our stiff fluid as a component with energy 
density ps, given by 

Ps = Psw{a/a 1Q y 6 . (7) 

where psw and aw are the stiff fluid density and scale 
factor, respectively, at T = 10 MeV (well before e + e~ 
annihilation). We use the Kawano [7l| version of the 
Wagoner [72|, [73[ nucleosynthesis code to derive the ele- 
ment abundances as a function of psw and of the baryon- 
photon ratio, r\. 

In Fig. 1 we compare the change in the element abun- 
dances produced by a stiff fluid with that from one extra 
neutrino species, for r\ in the range 5 — 7 x 10~~ 10 . The 
blue (solid) curve gives the standard BBN model with 
no additional energy density. The black (dashed) curve 
denotes the element abundances due to one additional 
two-component neutrino, while the red (dotted) curve 
gives the abundances due to a stiff fluid for which psw 
is chosen to produce the same effect on Y p as one extra 
neutrino (as can be seen in the top panel). 

As expected, when the stiff fluid density is adjusted to 
give the same effect on the 4 He abundance as an extra 
neutrino, the stiff fluid produces a much smaller effect 
on the deuterium and 7 Li abundances. This is because, 
as we have noted, the latter two element abundances are 
sensitive to the expansion rate at a much lower tempera- 
ture than is the 4 He abundance, and the stiff fluid density 



decays with expansion rate much more rapidly than does 
the contribution of an extra neutrino. 

The one minor surprise is that the stiff fluid increases 
the 7 Li abundance, while an additional neutrino de- 
creases it. This can be understood in terms of the re- 
sponse function in Ref. [70l |. The 7 Li response func- 
tion has a sharp trough just below T — 0.1 MeV, but 
it also has a shallow positive plateau for T > 0.1 MeV. 
The physical reason for this behavior comes from the fact 
that the dominant reaction for 7 Li production is 4 He + 
3 Hc — s> 7 Bc + 7 [67]]. An increase in the expansion rate 
at high temperatures (T > 0.1 MeV) gives a larger neu- 
tron abundance prior to nuclear fusion, enhancing the 
abundance of both 4 He and 3 He to make more 7 Be. An 
increase in the expansion rate just below T = 0.1 MeV, 
when nuclear fusion is occurring, gives less time for the 
fusion reactions to build into heavier elements; thus, the 
decrease in the 7 Li abundance in this case is accompa- 
nied by increases in the abundances of deuterium and 
helium-3. The stiff fluid samples the positive plateau at 
high T more strongly than the trough at low T, while the 
reverse is true for an extra neutrino. 

Once r\ is fixed by the CMB, the best constraint on 
the stiff fluid can be obtained from the 4 He abundance. 
WMAP7 gives r] = 6.2 x lO" 10 . For this value of 
r], we plot, in Fig. 2, the change in the primordial He 
abundance, AY p as a function of psw/ Prw> where psio 
is the density of the stiff fluid at T = 10 MeV (as in 
equation [7]), while prw is the standard energy density in 
relativistic particles at T = 10 MeV. As expected, we find 
that AYp is well- fit by a linear dependence on psw/ Prw, 
namely 

Ay p = 0.00024(p S io/Pflio)- (8) 

We find further that Eq. © is an excellent approxi- 
mation when rj lies in the range between 5 x 10~ 10 and 
7 x 10~ 10 . Of course, higher-order corrections produce 
a slight deviation from exactly linear behavior, as is ev- 
ident in Fig. 2 (see, e.g., Appendix F of Ref. 67] for a 
discussion of such corrections for AN V ). 

III. DISCUSSION 

The strongest constraints on a stiff fluid clearly come 
from the primordial 4 He abundance. This abundance re- 
mains at present somewhat uncertain (for recent discus- 
sions, see, e.g., Ref s. [75l - l77| ). Recent analyses by Izotov 
and Thuan [?5[ and by Aver, Olive, and Skillman [76| are 
both consistent with a central value of Y p = 0.256. Using 
this limit with the WMAP7 value of 77 = 6.2 x 10~ 10 , we 
obtain the bound 

Psw/ PRW < 30. (9) 

This result, however, should not be considered the 
main result of our paper, as the estimates of Y p are likely 
to improve with time. Rather, our main conclusion is 
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FIG. 1: Changes in primordial element abundances caused 
by a stiff fluid and by an extra neutrino, as a function of 
?7io = r\ x 10 , where r\ is the baryon- photon ratio. The 
black (dashed) curve denotes the element abundances due to 
one extra two-component neutrino, the red (dotted) curve 
denotes the abundances for a stiff fluid, and the blue (solid) 
curve denotes the unmodified standard case of 3 neutrinos and 
no stiff fluid. The density of the stiff fluid has been adjusted 
to produce the same effect on Y p as one extra neutrino, as can 
be seen from the top panel. 



embodied in Eq. ((5J, a result that can be used to pro- 
vide an upper bound on the stiff fluid density for any 
estimate of the primordial 4 He abundance. 

We can exploit the fact that AY p depends linearly on 
both AiV„ and psio/pmo to derive a "pivot tempera- 
ture," T*, at which equal contributions from relativistic 
energy density or from a stiff fluid will produce equal 
changes in 4 He. In other words, suppose that we have a 
particular value of AY p . This will correspond to a partic- 
ular value of AN U in Eq. §|, or to a particular value of 
Psio/pmo in Eq. ©. In this case, the additional energy 
densities in radiation or in the stiff fluid will be equal at 
a single temperature T»: 

A PR (T*) = p s (T*). (10) 

By comparing Eqs. © and ([5]), we find that 

T* = 0.55 MeV. (11) 

We emphasize that this is a purely heuristic result. It 
is not based on the assumption that the helium abun- 
dance depends only on the expansion rate at a single 



temperature; we have seen that it does not. Indeed, it 
is not surprising that T* lies in between the two peaks 
in the 4 He response function. Note further that T* cor- 
responds, strictly speaking, to the neutrino temperature, 
rather than the photon temperature. The neutrino tem- 
perature scales exactly as 1/a, while the photon tem- 
perature, even near 0.55 MeV, has experienced a small 
increase due to e + e~ annihilation. Since most bounds 
from BBN on additional energy density are expressed in 
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FIG. 2: The change in the primordial 4 He abundance, AY P , 
produced by a stiff fluid with density ps = psio(«/iio) _6 , 
where psio is the stiff fluid density at a temperature of 10 
MeV, and aio is the scale factor at this temperature. The 
horizontal axis shows psio/pmo, where pmo is the total en- 
ergy density in relativistic particles (i.e., in the standard cos- 
mological model without a stiff fluid) at 10 MeV. Blue (solid) 
curve is numerical result, and red (dotted) curve is the fit 
from Eq. ((HJ. These results are for a baryon-photon ratio of 
77 = 6.2 x 10~ 10 . 



terms of limits on additional radiation density (or equiv- 
alently, additional neutrinos) from the primordial helium 
abundance, our result for the pivot temperature can be 
used to convert these bounds into limits on an additional 
stiff fluid component. 



IV. ACKNOWLEDGMENTS 

We thank E. Linder, G. Steigman and and R. Stiele for 
helpful comments on the manuscript. We thank E. Sari- 
dakis for useful discussions. R.J.S. was supported in part 
by the Department of Energy (DE-FG05-85ER40226). 



[1] T.P. Walker, G. Steigman, D.N. Schramm, K.A. Olive, 
and H.S. Kang, Ast rophys. J. 376, 5 1 (1991). 

[2] E. Komatsu, et al. arXiv:1001.4538 ; D. Larson, et al. 
[arXiv: 1001.4635] . 

[3] M. Hicken et al, Astrophys. J. 700, 1097 (2009) 



|arXiv:0901. 48041 [astro-ph.CO]] . 
[4] W. J. Perci val et al, Mon. N ot. Roy. Astron. Soc. 401, 

2148 (2010) arXiv:0907.1660l [astro-ph. CO]]. 
[5] L. Samushia, G. Chen and B. Ratra, larXiviOTOe.Tglal 

[astro-ph] . 



5 



1347 



645 



S. Ettori et al, arXiv:0904.2740 [astro-ph.CO]. 
Y. Wang, Phys. Rev. D 78, 123532 (2008). 
L. Samushia and B. Ratra, Astrophys. J. 714, 
(2010), [arXiv:0905.3836] . 

A. Refregier, Ann. Rev. Astron. Astrophys. 41, 
(2003) [arXiv:astro-ph/0307212] . 
J. A. Tyson, G. P. Kochanski and I. P. D ell' Anto nio, A s- 
trophys. J. 498, L107 (1998) [arXiv:astro-ph/9801193] . 
J. M. Benko, G. A. Bakos and J. Nuspl, Mon. Not. Roy . 
Astron. Soc. 372, 1657 (2006) [arXiv:astro-ph /0608409] . 
S. W. Allen, A. C. Fabian, R. W. Schmidt and H. Ebel- 
ing, Mon. Not. Roy. A stron. Soc. 342, 287 (2003) 
[arXiv:astro-ph/0208394] . 
F. Zwicky, Helv. Phys. Acta 6, 110 (1933). 
V. C. Rubin and W. K. J. Ford, Astrophys. J. 159, 379 
(1970). 

V. C. Rubin, N. Thonnard and W. K. . Ford, Astrophys. 
J. 238, 471 (1980). 

E.J. Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. 
Phys. D 15, 1753 (2006). 
J. L. Feng. larXiv: 1003.09041 [astro-ph.CO], 
S. Dutta, S. D. H. Hsu, D. Reeb and R. J. Scherrer, 
Phys. Rev. D 79, 103504 (2009) |arXiv:0902. 4699 [astro- 
ph.CO]]. 

Zeldovich, Ya. B. Sov. Phys. JETP 14, 11437 (1962). 
M. Joyce, Phys. Rev. D 55, 1875 (1997) 



arXiv:hep-ph/9606223 



M. Joyce and T. Prokopec, Phys. Rev. D 57, 6022 (1998) 
arXiv:hep-ph/9709320 . 



M. Kamionkowski and M.S 
3310 (1990). 

P. Salati, Phys. L ett. 
|arXiv:astro-ph/0207396] . 
C. Pallis, Nucl. Phys. B 751 



B 571, 



Turner, Phys. Rev. D 42, 
(2003) 
(2006) 

and J. Rodriguez- 



121 



129 



|arXiv:08124287l 
084008 (2009) 



arXiv:hep-ph/0510234 . 
M.E. Gomez, S. Lola, C. Pallis, 
Quintero, JCAP 0901, 027 (2009). 

C. Pallis, Nucl. Phys. B 831, 217 (2010) |arXiv:0909.3026l 
[hep-ph]]. 

D. J. H. Chung, L. L. Ever ett and K. T. Ma tchev, Phys. 
Rev. D 76, 103530 (2007) |arXiv:0704.3"285l [hep-ph]]. 
R. Stiele, T. Boeckel and J. Schaflher-Bielich, Phys. Rev. 
D 81, 1235 13 (2010), |arXiv:1003.2304| . 

P. Horava, larXiv:0811.2217| . 
P. Horava, JHEP 0903, 020 (2009) 
[hep-th]]. 

P. Horava, P hys. Rev. D 79, 
|arXiv:0901.3775l [hep-th]]. 
P. Hofava, arXiv:0902.3657 [hep-th]. 
G. Calcagni, larXiv:09 04.0829 [hep- th]. 

E. Kiritsis and G. Kofinas. la7Xiv:0904. 13341 [hep-th] . 

T. P. Sotiriou, M. Visser and S. Weinfurtner, JHEP 
0910, 033 (2009) |arXiv:0905. 27981 [hep-th]]. 
C. Bogdanos and E. N. Saridakis, larXiv:0907.To36l [hep- 
th]. 

S. Carloni, E. Elizalde and P. J. Silva, I arXiv: 0909. 22191 
[hep-th]. 

G. Leon and E. N. Saridakis, JCAP 0911, 006 (2009) 
|arXiv:0909.357T1 [hep-th]] . 
S. Mukohyama, K. Nakayama, 
and S. Yokoyama, Phys. Lett. B 
|arXiv:0905.0055l [hep-th]] . 
S. Mukohyama, JCAP 0906, 
|arXiv:0904.2T90l [hep-th]] . 



F. 

679, 



Takahashi 
6 (2009) 



001 (2009) 



[41 
[42 
[43 

[44 
[45 
[46 
[47 



[49 
[50 

[51 

[52 
[53 
[54 

[55 
[56 
[57 

[58 
[59 

[60 

[61 

[62 

[63 

[64 

[65 

[66 

[67 
[68 

[69 

[70 

[71 
[72 

[73 
[74 

[75; 

[76 
[77 



R. Brandenberger, Phys. Rev. D 80, 043516 (2009) 
|arXiv:0904.2835l [hep-th]] . 

Y. F. Cai and E. N. Saridakis, JCAP 0910, 020 (2009) 
|arXiv:0906. 17891 [hep-th] ] . 

A. Wang and Y. Wu, JCAP 0907, 012 



B 684, 194 



(2009) 
(2010) 



1005, 013 (2010) 
1001, 013 (2010) 



015 (2009) 



arXiv:0905.4117 [hep-th]]. 
M. Minamitsuji, Phys. Lett 
[arXiv:0905.3"892l [astro-ph.CO]]. 

A. Ali, S. Dutta, E. N. Saridakis and A. A. Sen 
larXiv:1004.2474l [astro-ph.CO]. 
S. Dutta and E. N. Saridakis, JCAP 
[arXiv:1002. 33731 [hep-th]] . 
S. Dutta and E . N. Saridakis, JCAP 
[arXiv:0911. 14351 [hep-th]] . 

D. Bias , O. Pujolas and S. Sibiryakov, JHEP 0910, 029 

(2009) _arXiv:0906.3046 [hep-th]]. 

A. Kobakhidze, arXiv:09 06340T1 [hep-th] . 
M. Li and Y. Pang, JHEP 0908 
[arXiv:0905.275T1 [hep-th]] . 
M. Henneaux, A. Kleinschmidt and G. L. Gomez, Phys. 
Rev. D 81, 064002 (2010) |arXiv:0912.0 399 [hep-th]]. 
J. M. Pons and P. Talavera, arXiv: 1003.3811 [gr-qc]. 
K. Koyama and F. Arroia. larXiv :0910.1998 [hep-th]. 
A. Papazoglou and T. P. Sotiriou, larXiv:0911.T299l [hep- 
th]. 

I. Kimpton and A. Padilla. lajXivT l003.5666 [hep-th]. 
J. Bellorin and A. Restuccia, arXiv: 1004.0055 [hep-th]. 
L. Fernandez-Jambrina and L. M. Gonzalez-Romero, 
Phys. Rev. D 66, 024027 (200 2) [arXiv:gr-qc/0402"TT9] . 
N. Dadhich, [arXiv:gr-qc/9406034| 

M. Mars, Phys. Rev. D 51, 3989 (1995) 
[arXiv:gr-qc/0201078j . 

L. Fernandez-Jambrina, Class. Quant. Grav. 14, 3407 
(1997) [arXiv:gr-qc/0404017| . 

F. Donato, N. Fornengo, and M. Schelke, JCAP 0703, 
021 (2007) [arXiv:hep-ph/0612374| . 

M. Drees, H. Iminniyaz, and M. Kakizaki, Phys. Rev. D 

76, 103524 (2007) [arXiv: 0704. 1590] , 

A. Arbey and F. M abmoudi, Phys. Lett. B 669, 46 (2008) 

[arXiv:0803.0741| . 

A. Arbey and F. Mahmoudi, JHEP 1005, 051 (2010). 
[arXiv:0906.0368| 

E. Masso and R. Rota, P hys. Rev. D 68, 123504 (2003). 
[arXiv:astro-ph/0302554] 

S. Carroll and M. Kaplinghat, Phys. Rev. D 65, 063507 
(2002). 

PD. Serpico, et al, JCAP 0412, 010 (2004). 

G. Steigman, An n. Rev. Nucl. Part. Sci. 57, 463 (2007). 
[arXiv:0712.1100| 

F. Iocco, G. Mangano, G. Miele, O. Pisanti, and P.D. 
Serpico, Phys. Rept. 472, 1 (2009) [arXiv: 0809.0631] . 
C. Bambi, M. Giannotti, and F.L. Villante, Phys. Rev. 
D 71, 123524 (2005), [arXiv:astro-ph/0503502] . 

L. Kawano, Fermilab-pub-92/04-A (1992). 

R.V. Wagoner, W.A. Fowler, and F. Hoyle, Astrophys. 

J. 148, 3 (1967). 

R.V. Wagoner, Astrophys. J. 179, 343 (1973). 

J.P. Kneller and G. Steigman, New J. Phys. 6, 117 

(2004). 

Y.I. Izotov and T.X. Thuan, Astrophys. J. Lett. 710, L67 

(2010) . 

E. Aver, K.A. Olive, and E.D. Skillman, 
[arXiv:1001.52T8] . 

G. Steigman, JCAP 1004, 029 (2010), [arXiv:1002.3604] . 



